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Abstract 

We propose a new Complex Diffusion Monte Carlo (CDMC) method 
for the simulation of quantum systems with complex wave function. In 
CDMC the modulus and phase of wave function are simulated both in 
contrast to other methods. We successfully test CDMC by the simulation 
of the ground state for 2D electron in magnetic field and 2D fermions- 
anyons in parabolic well. 

1 Introduction 

In the quantum mechanics, there are several cases where the wave function 
is essentially complex(has imaginary part). The examples are 2D system of 
electrons in the external uniform magnetic field when vector potential has a 
central-symmetrical form and system of anyons. The simulation of these sys- 
tems is impossible by using well-known Green Function Monte Carlo (GFMC) 
method (see review [Q]) because GFMC requires the reality of wave function 
which is considered as the probability weight during stochastic process. 

There have been several attempts [||, ||, |j to construct a Monte Carlo 
method for the simulation of quantum systems with complex wave function. 
Authors of Ref . § tried to apply GFMC to this problem and met essential diffi- 
culties. Another method Q has been developed recently by using an algorithm 
without branching that gave the increasing of the systematical statistical errors 
as a function of the simulation time. 

The basic difficulty of numerical simulations of fermions is essentially the 
same as for the systems with complex wave function. Their wave function 
can change the sign and therefore, can not be used also as the probability 
weight. For the simulation of fermionic systems a widely used fixed node Monte 
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Carlo method has been developed ^ (see also ||). Recently the constrained 
path Monte Carlo method was proposed J?]]. In both methods one assumes the 
restriction on the random walks connected with the uncertainty in the space 
localization of the wave function node surface. CDMC method is very close 
to the fixed phase diffusion Monte Carlo method Q which was applied to 2D 
electrons in magnetic field. In the framework of diffusion Monte Carlo method 
one simulates only modulus of the system wave function. The phase of this 
wave function is not simulated but treated to be fixed and equal to the phase of 
Laughlin's wave function ||. Here we propose a new CDMC method including 
also the simulation of the phase of the system wave function |l0| . In the 2D 
space we have unique tool - so-called anyons, particles (bosons or fermions) 
with additional gauge interaction which provide required statistical properties 
|Tl| l . So, by tuning the coupling constant we may arrive to the fermions starting 
from the bosons. It looks very attractive to apply the anyons for the simulation 
of 2D fermions because statistical vector potential that gives anyonic property 



is relatively smooth and continuous |Ll], 12 1. The remaining main problem is 
that the wave function of anyonic system contains an imaginary part. 

The first test of CDMC reproduces the ground state energy and its degen- 
eracy on the orbital quantum number m of electron in magnetic field. Also the 
simulated wave function is in exact correspondence with analytical prediction 
(Fig.(jlj)). As the second test of CDMC, we simulate ground state energy for 
2D fermions and anyons in parabolic well. The system of 2D fermions has a 
well-known analytical solution for this energy as a function of the number of 
fermions. The simulation is done for the fermionic systems with the number of 
particles from two to ten. We have found a good agreement of the numerical 
results with analytical ones for the numbers of the particles great than four 
(Tables (jl|), @). This observation of simulation allows us to hope that the 
CDMC is a good tool for the simulation of 2D fermion systems with big num- 
ber of particles. The simulated ground state energy for anyons (Table (|^)) we 
compare with variational calculation 13]. 



2 Ground state of 2D electron in magnetic field 

Hamiltonian of the electron with A = [Hr\/2 (in the gauge pA—Ap = divA = 0) 
has a form: 

Here p = —ihV, V = (d/dx, d/dy). 

Let us introduce a complex distribution function 

f(r,t) = y* T (r)V(r,t). (2) 
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Here ^^(r) is a complex conjugated trial wave function of electron in the mag- 
netic field. Schrodinger equation for VP(f, t) with imaginary time ( expressed in 
H units ) is: 

-^ = (ff-W,t). (3) 
The equation for the distribution function / = f(f, t) is: 



df 
dt 



Here we have 



-DAf + DV(f ReF Q {r)) + i[V(Df ImF Q {r)) - 
MAVf] + (E L (r)-E T )f. (4) 



F Q (r)=2** T 1 (r)V*Ur), (5) 
£ L (f) = ^-^fO^'^CrO, (6) 

H' = — (p - ^i") 2 (7) 
ft 2 

where we putted p = —ih"V, D = , A = V 2 and took into account 

v 2M' 

F Q (r ) = ReF Q (r) + ilmFq (r ) (8) 

when ^y(r) is complex. 

Following || , in the limit r — > for the time step of equation (Q) we assume 
FQ^f) = Fq^') where f corresponds to time point t + r and r ' to t. 

Let us introduce a new quantity 

A Q (f, r ') = \lmF Q (r ') - ^f/i?). (9) 
Then at r — ► the Green function of equation (Q) has a form: 
G(r,r ';r) = Gi(f,f ';r)exp [-r(E L (f) - E T )\ x 



x exp 

where 



il Q (r,r ')(r-r ' - DrReF Q (f ')) (10) 



, exp[Dr^(r,r')] 
Gl(r ' r ;T) = 4^7 eXP 



(f - f ' - DTReF Q (f')) 2 



ADt 



(11) 
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The Green function G(f,f ';r) (10) and distribution function / (|2|) are the 
complex functions and related by usual integral equation: 

f(f,t + r) = J df'G(f,r';r)f(f f ,t). (12) 

From ( |i~2| ) the modulus and phase of / at t + r are determined by the modulus 
and phase of Green function and of ones of / at t. 

If ^y(r) is complex then the energy E^r) is a complex too. Therefore, the 
last two exponents in ([TOD have a form: 



exp {-T{ReE L (r) - E T )\ 



x exp 



iA Q {r,r')(r-r'-DTReF Q {r'))-iTlmE L (r) . (13) 



CDMC algorithm |] is: 

1 . We prepare iV c initial configurations r in which the electron has the random 
and uniform distribution in a position. 

2. We perform the quantum drift and diffusion of the particle: r*/% = ' + 
DrReFQ^k ') + x, for example, from k -th configuration. Here x is Gaussian 
random number having zero mean value and the dispersion 1\J Dr. 

3. The transition into new space point is accepted with probability P(f ' — > 

r,r) ^ min(l, W(r,r ')) where W(r,r>) = ^^^7^ - <h<?f V) 

is given by (|Tl|). If the transition of electron is accepted then in accordance 
with (|D|) it has a new phase; if not then the electron keeps its old phase. 

4. After changing the electron position from /^-configuration into new space 
point, one calculate ReEiirk), ImEL(rk) and other quantities for measure- 
ments. 

5. By using the first exponential factor in ( |l3"| ) we calculate the multiplicity Mk 
according to the formula Mk = exp[—T(ReEL(fk) — Et)\- If is not integer 
we add an uniformly distributed random number between and 1 to it and 
take Mk equal to nearest integer. 

6. If Mk ^ then Mk copies of new A:-th configuration are put in the list 
of the new ./V configurations which becomes the initial one at the next step r 
in the integration of the diffusion equation. If Mk = then there is no fc-th 
configuration in the list of new N configurations for the next time step r. All 
quantities of interest such as ReEi(fk), IraE^ifk) etc. are multiplied by the 
factor Mk when we calculate their mean values. 

7. We repeat the steps (2) - (6) until all N c configurations are not overlooked 
and electrons on these configurations are not simulated on the displacement. 

8. We calculate a mean energy and other mean quantities on the N configu- 
rations obtained at the step (6) at the time step r in accordance with formula 
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(!)• 

9. We repeat the points (1) - (8) of algorithm an integer number of time steps r 
of integration of diffusion equation. After that one determine the mean values of 
quantities ReE, ImE and others on this integer number too. New value of 

is redetermined in accordance with the formula (Ex)new = [{Ex) old + ReE]/ 2 
with the assumption ImE << ReE. The integer number of time steps r rep- 
resents one time block At. 

10. Every time block At has the N c initial configurations. The list of these 
configurations is filled by randomly chosen configurations from just ended time 
block. The random choice of configurations consists from two steps: 

a) a random choice of configuration in the every time step r in the time block 
At; 

b) in the case N c larger than number of the time steps r in the time block At, 
other configurations are taken randomly from the last time step r. 

In this manner filled list of the N c configurations will be initial one for the next 
time block At. 

11. The big number repeating of time blocks At essentially decreases a cor- 
relation between configurations in neighbor time blocks and provides a right 
calculation of mean quantities. 

The calculation of the mean quantities by complex distribution function 
f(R,t) at the time t of the running process is occurred by the formula Q: 

M 

J2^MHRi(t)W(Ri(t)) 

<F(t)> = *=L_ (14) 

exp[ia(Ri(t)] 

i=i 

M 
i=l 



M 

y-ya G (-Ri(t),£i(*-r)) 
i=l 



(15) 



Here Ri(t) is fi, r*2, nv> N is the number of particles and M is the number 
of configurations at time point t. The particle coordinates ri,r*2, ...,rjv in (111), 
(|l^) are weighted with probability \f(Ri,t)\ and a(Ri(t)) is a phase of the 



distribution function /(i?j,t)and ac is the phase of Green function (13). 

The step from Eq. ( |l4l) to Eq. (15) is rather nontrivial and can be 
proved with account the properties of the phase of Green function (|l3|) oiq. 
Because at r — > it has an asymptotic ac ~ t 1 / 2 (see ||) and oscillates ran- 
domly around zero. In the following, the comparison of the results of numerical 
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Figure 1: Left panel: The real part of the ground state energy ReE as function of 
the number T of time blocks At for different orbital quantum numbers m. Solid line, 
marked solid line and dashed line correspond m = 0, m = 4 and m = 8, respectively. 
Here a = (3 — 7 = 1.0. Right Panel: The simulation result for ground state spatial 
distribution of the electron in uniform magnetic field at m — 13. The length unity is 
a c , initial number of points N c — 500 and parameters a = j3 = 7 = 1. 



calculation of ground state energy of 2D fermions by the formula (O) (Table 
(||)) and by the formula (|l^) (Table (|l])) confirms this suggestion. 

The trial wave function for the ground state of 2D electron in magnetic field 

is: 

* {f) = — exp -7 — . (16) 

a c \ a c J y 4a c J 

C - normalizing constant, a c = (h/Mw c ) l l 2 - magnetic length and w c = 
\e\H/Mc - cyclotron frequency, x and y - components of electron coordinate, 
m is the orbital quantum number, and a, (3 and 7 are arbitrary numerical con- 
stants. When a = (3 = 7 = 1 ( |l6| ) is an exact ground state wave function. We 
take a c , hw c /2, and 2/hw c as length, energy and time units, respectively. A 
number of time blocks At is ten. For every At the initial number of configura- 
tions N c = 1000 and number of time steps r is equal 100. Everywhere r was 
equal 0.01. 

Left panel of Fig. (JTJ) presents ReE as functions of number time blocks T. 
Right panel of Fig.((l|) is the ground state spatial distribution of electron after 
simulation for m = 13. From this Fig. we see good correspondence with 



exact solution [15] for the spatial distribution which is a ring with mean radius 
(2m + l) 1 / 2 and wide 1 in a c magnetic length units. 
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3 Ground state for 2D fermions and anyons in parabolic well 

Hamiltonian of anyons in the 2D parabolic well has a form: 

1 N A/ 2 W 

h = ^ EG* + ^)) 2 + ^ E *?■ 

1=1 1=1 

Here M is the mass of particle, pi = —ihVi where V, = (d/dxi,d/dyi), uj Q is a 
characteristic frequency of free particles in parabolic well and f{ is radius vector 
of i-th particle, N is number of particles. 

Vector potential for anyons A(fi) jll], 12] in (|l7|) is 

N -> -> 



j>i 



Here z is the unit vector perpendicular to 2D plane and v - anyonic fractional 
statistics factor (spin of the anyon). In the bosonic representation of anyons 
v = gives noninter acting bosons and v = 1 gives free fermions. 
We introduce a complex distribution function 

f(R,t) = y* T (R)y(R,t) (19) 

where R stands for the coordinates of all particles and the wave function ^f(R, t) 
satisfies a Schrodinger equation (|3]) with hamiltonian fll7|). 

For a bosonic representation of anyons we take (conjugated) trial wave func- 
tion in the form: 

JV 

^(i?) = n*r(n). (20) 

8=1 

The distribution function / = f(R,t) satisfies a diffusion equation: 
Qf N N 

-^ = - D Y, A *f + D Y, Vi(/ReF Q (r-)) + 

i=l i=l 
N j- 
+i^ l (DfImF Q (f l )) - —A(n)V t f] + (E L {R) - E T )f. (21) 

i=l 

Here D = h 2 /2M and Aj = Vf . When a time step integration of diffusion 
equation ( |2l| ) r goes to zero, the Green function for this equation has a form: 

N 

expprj^for* ')] 

G(R,R';r) = x 
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x exp 



N 

£( 

i=l 



DrReFg (Fi ')) 2 



4L>r 



cxp 



-t{E l {R)-E t ) 



x exp 



A Q (n,fi ')(F - F ' - DrReFg^ ')) 



(22) 



In the expressions (^) and ( p2| ) .Fq(F) determines by formula (a) where radius 
vector and operator nabla have index i, the energy El(R) determines by formula 
(||) with wave function ^^(R) and hamiltonian H' (the hamiltonian ([l?]) with 
changing of sign in circle brackets (see (0) for determination of H')). i*Q(F) for 
the complex wave function ^^(F) has the real and imaginary parts. 
In (|22|) we have introduced a new quantity 



^Q(F,F ') = ^ImF Q (F ') - ^A(F). 



(23) 



In the expressions (^) and (^) vectors i? and F correspond to the time point 
t + r and vectors R ' , f ' to the time point £. 

As the energy Ei{R) is complex, so the expression for two last exponents 
in (22) has a form: 



exp 



-T(ReE L (R) - Et) 



x exp 



iV 



: ]T A Q (F, F ')(F - F ' - DrReF Q (F ')) " iTlmE L (R) 



i=i 



(24) 



mean magnetic field [16] with vector potential 



The ground state trial wave function for one anyon in parabolic well and in the 

(25) 



1 -* 



A(r) = pirhv(zx r) = -B x F 
which is generated by the average density of particles p can be chosen as: 



^r(F) = Cexp -a 



(A + y]) 

2Rl 



exp 



Aa 2 H 



(26) 



Here B = lixphvz is fictitious uniform mean magnetic field. Constant C is 
for normalization and R Q = (Ti/Moj ) 1 I 2 . Density p determines by formula 
p = l/wr 2 , where r D is the mean distance between particles. We assume (see 
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Figure 2: Square of module of one particle ground state wave function for 6 fermions 
^(r)! 2 = J \^>(r — fi, T2, ■ ■ ■ , fe)\ 2 df^df's ■ ■ • dr§. Left panel: result of the simulation. 
Right Panel: calculations with Slater determinant. 



llOfl ) that r Q is equal to R a . We express the energies in terms of hu> and 
the lengths in terms of r Q . We calculate anyons ground state energy at v = 
0.2, 0.4, 0.6, 0.8. The total simulation time T is equal 40 time blocks At. Anyons 
become fermions at v = 1. 

For every At N c = 1000 and the number of time steps r is 30. 

We take here optimal time step To in accordance with |1C]. We find that 
the imaginary parts of the simulated mean quantities for all numbers fermions 
and anyons are essentially less than their real parts. 

The results of the simulation for 2D fermions-anyons are given in tables 
(|) - ©. In tables (0) and (|) exact energies E were calculated by occupying 
the available lowest states of 2D parabolic well while the minimization of the 
simulated energy over parameters a gives the E$im- The comparison of the 
square of module of simulated one particle ground state wave function for 6 
fermions in parabolic well with the square of the same module calculated by 
Slater determinant is given in the figures @. From tables (|l]) and (]2|) we 
see that for four and bigger numbers of 2D fermions we have good simulation 
results. Good results for the simulation of ground state for the electron in 
magnetic field(see figures (|l|)) and for the systems of 2D fermions in parabolic 
well (tables (||), (||), figures (§)) (for the number of particles bigger than four) 
allows us to hope that CDMC provides the successful simulation of 2D fermionic 
systems with Coulomb interaction. 

We are grateful to S.Fantoni, M.-P. Lombardo, G.Ortiz and G.Senatore 
for useful discussions. One of us (B.A.) would like to thank organizers for the 
support of his participation in Quantum Monte Carlo 2001 workshop in Trento. 



10 



B. Abdullaev, M. Musakhanov, A. Nakamura 



N 


E 


Esim 


r 




a 


To 


2 


3 


1.842 
±1.297 • 10" 3 


0.935 
±6.994 • 10" 4 


0.8852 
±1.319 • 10" 3 


-0.593 


0.08 


4 


8 


7.399 
±4.092 • 10" 2 


1.236 
±1.844- 10" 2 


1.919 
±4.562 • 10" 2 


- 0.3 


0.04 


6 


14 


14.492 
±3.136 • 10" 2 


1.619 
±4.261 • 10" 3 


3.275 
±1.564- 10" 2 


-0.092 


0.02 


8 


22 


22.383 
±8.569 • 10" 2 


2.020 
±5.821 • 10" 3 


5.112 
±2.790 • 10" 2 


-0.046 


0.015 


10 


30 


30.905 
±2.010 • 10" 1 


2.601 
±2.839 • 10" 2 


8.599 
±2.095 • 10" 1 


-0.084 


0.005 



Table 1: The ground state energies for the fermions in 2D parabolic well. Here N - 
numbers of particles, E - analytically calculated ground state energies, Esim - results 
of the simulation calculated by formula ( |l5| ) (both in hu units), r - the simulated 
mean radius (in r a units), r 2 - the simulated mean square radius (in r 2 , units), a - the 
numerical parameters in the wave function ( p6| ) that give a minimum Esim and tq - 
the optimal time steps. All simulated quantities and their deviations from mean values 
were averaged over 30 last time blocks At when Esim and mean population number 
N p have had relatively stable values (N p as 1000). 
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N 


E 


Esim 


r 




a 


TO 


2 


3 


1.842 
±1.192 ■ 10" 3 


0.935 
±6.974 • 10" 4 


0.8847 
±1.418 • 10" 3 


-0.593 


0.08 


4 


8 


7.402 
±2.374 • 10" 2 


1.242 
±1.202 • 10" 2 


1.936 
±3.149 • 10" 2 


- 0.3 


0.04 


6 


14 


14.493 
±4.160 • 10" 2 


1.612 
±5.206 • lO" 3 


3.245 
±2.219 • 10" 2 


-0.092 


0.02 


8 


22 


22.358 
±6.639 • 10" 2 


2.011 
±6.379 ■ 10" 3 


5.075 
±3.065 • 10" 2 


-0.046 


0.015 


10 


30 


30.896 
±2.303 • 10" 1 


2.598 
±4.129 • 10" 2 


8.574 
±2.941 • 10" 1 


-0.084 


0.005 



Table 2: The ground state energies for the fermions in 2D parabolic well. Here Esim 
- results of the simulation calculated by formula (|l4|), the determinations of other 
quantities are the same as in table (m) and N p w 1000. 



V 


E v 


Esim 


r 




a 


TO 


0.2 


32.00 


30.815 
±3.222 


1.353 
±7.078 • 10" 3 


2.282 
±2.430 • 10" 2 


0.500 


0.001 


0.4 


42.33 


40.525 
±0.268 


2.081 
±3.611 • 10" 2 


5.508 
±2.040 • 10" 1 


0.500 


0.001 


0.6 


50.59 


48.371 
±1.605 


2.260 
±5.124- 10" 2 


6.185 
±2.771 • 10" 1 


0.500 


0.0075 


0.8 


57.68 


54.119 
±1.142 


3.085 
±7.307 • 10" 2 


11.984 
±6.544 • 10" 1 


0.177 


0.005 


1.0 


64.00 


65.401 
±0.914 


3.646 
±6.379 • 10" 2 


16.872 
±7.028 • 10" 1 


-0.050 


0.005 



Table 3: The ground state energy as function of anyonic factor v for 16 anyons in ID 
parabolic well. Here v - anyonic factor, E v - a variational ground state energies ||l3f 
(in hue, units) and N p w 1000 (except N p 2000 at v = 0.6). The determinations of 
other quantities are the same as in tables (Eh and (§). 



